Large Scale Structure as a Probe of Gravitational Slip 
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A new time-dependent, scale-independent parameter, w, is employed in a phenomenological model 
of the deviation from General Relativity in which the Newtonian and longitudinal gravitational po- 
tentials slip apart on cosmological scales as dark energy, assumed to be arising from a new theory 
of gravitation, appears to dominate the universe. A comparison is presented between w and other 
parameterized post-Friedmannian models in the literature. The effect of w on the cosmic microwave 
background anisotropy spectrum, the growth of large scale structure, the galaxy weak-lensing cor- 
relation function, and cross-correlations of cosmic microwave background anisotropy with galaxy 
clustering are illustrated. Cosmological models with conventional maximum likelihood parameters 
are shown to find agreement with a narrow range of gravitational slip. 



I. INTRODUCTION 

The quest for the source of the cosmic acceleration has 
lead to speculation that the proper theory of gravitation 
departs from General Relativity (GR) on cosmological 
scales. That is, gravitation may be well described using 
Einstein's theory within the solar system and the envi- 
rons of the galaxy, but a different theory is required on 
the scale of the Hubble length (e.g. Ref. [l[). There are 
numerous examples of a theory capable of producing a 
late-time acceleration, all of which introduce new gravi- 
tational degrees of freedom, with wide-ranging implica- 
tions for observable phenomena 0, Q. Given this pos- 
sible abundance in new physics, it is important to come 
up with tests that can distinguish between the effects 
of dark energy and those of modified gravity. Though 
late-time accelerated cosmic expansion is the principal 
indicator that a new "dark" physics is needed, it is not 
the only test such physics must satisfy. A successful cos- 
mology must also agree with measurements related to 
the growth of perturbations as probed by the cosmic mi- 
crowave background, clustering of large-scale structure, 
and weak lensing deflections of light. 

The concordance cosmology, ACDM within GR, gives 
very specific predictions for the cosmic expansion and 
the evolution of inhomogeneities. It has been proposed 
to evaluate alternative theories of gravitation by testing 
for violations of these predictions (e.g. Ref. [J]). Ef- 
forts along these lines consist of fitting separate ACDM 
parameters to tests of structure growth and cosmic ex- 
pansion or constraining a phenomenological descrip- 
tion of the linear growth rate of density perturbations 
with galaxy clustering surveys @, 0, H, Q • To within the 
stated uncertainties, these current results are all consis- 
tent with ACDM. The absence of any realistic theory of 
the cosmological constant, however, has prompted many 



to question whether gravity itself is the culprit. 

An exhaustive study of departures from GR as an 
explanation for dark energy phenomena is difficult, as 
Einstein's theory represents a mere island in a sea of 
possible gravitational theories. Yet, a common feature 
within a broad range of such theories is a decoupling 
of the perturbed Newtonian-gauge gravitational poten- 
tials and ip, defined by the perturbed Robertson- Walker 
line-element 



ds 2 = a 2 [- (1 + 2-0) dr 2 + (1 - 2</>) dx 2 



(1) 
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using the notation and convention of Ma & Bertschingcr 
[Tol | ■ T° gi ye a physical sense of these potentials, ip enters 
the Newtonian limit of the equation of motion, x = —Vip, 
and 4> enters the Poisson equation V 2 4> = ~AnGa 2 8p. 

Whereas GR predicts ip = <f> in the presence of non- 
relativistic matter, a gravitational slip, defined as ip ^ 
4>, generically occurs in modified gravity theories. This 
inequality means that the gravitational potential created 
by a galaxy cluster is not the same potential responsible 
for the geodesic motion of the constituent galaxies. For 
primordial cosmological perturbations, the potentials are 
not completely free, however, as there exists a constraint 
equation, valid under general assumptions in the long- 
wavelength limit [Tll |. Hence, a new relation between 
these potentials is a launching point for investigations of 
cosmological manifestations of modified gravitation. 

With a view towards testing modified gravitation 
against large-scale structure, the relation between 4> and 
■tp has been examined for scalar-tensor theories [h|, 
generalized gravitational Lagrangians or f(R) theories 
0, [3 , a tensor- vector-scalar model of gravity fill , the 
Dvali-Gabadze-Porrati (DGP) model [la, [TtI Il8l . and 
a Lorentz-invariance violating massive gravity [1 91 ] . Phe- 
nomenological relations between the potentials have been 
widely investigated [H HI, [H, H, IM H, [Hj| . A gen- 
eral set of post-Friedmannian parameters has also been 
proposed [271 ] , which reproduces the evolution of (j> and vb 
on both sub- and super-horizon scales for DGP and f(R) 
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gravities; the relation between these parameters and cos- 
mic acceleration has been discussed in Ref . [28[ . Much of 
the analysis has been formal, with an aim towards future 
tests. Few efforts have been made to place constraints 
using current data. 

In this paper we explore the parameterized post- 
Friedmannian (PPF) description of gravitation intro- 
duced in Ref. [29[ (hereafter CCM). CCM posit a mod- 
ified theory of gravitation that produces a ACDM- 
cquivalcnt background with perturbations such that ip = 
(1+vj)4>, where vo is taken to be scale-independent. Here, 
we correct an important error in CCM, and explore the 
impact on a wider range of cosmological phenomena. In 
Sec. [TT]we review the model for the time-dependent quan- 
tity vo and the description in terms of cosmological pa- 
rameters. We compare the PPF formalism to similar 
attempts at modifying GR in the literature. In Sec. IIIII 
we give the procedure for evolving cosmological perturba- 
tions and describe modifications to the CMBfast software 
[30l | ■ In Sec.[lV]we explore the influence of vo on the cos- 
mic microwave background (CMB) anisotropy, correct- 
ing several errors in the numerical calculations described 
in CCM. We show the effect ofw^Oon the Wilkin- 
son Microwave Anisotropy Probe (WMAP) [3l[ best-fit 
cosmology and its consistency with current data. Us- 
ing WMAP best-fit cosmological parameters, we simi- 
larly demonstrate the effect of vo ^ on perturbation 
growth in Sec. [Vj weak lensing in Sec. IVI1 and the in- 
tegrated Sachs- Wolfe effect in Sec. IVIII We attempt a 
synthesis of these results in Sec. lVIIll and present a final 
discussion in Sec. IIXI 



II. PARAMETERIZATION 

Non-Einstein gravitation gcncrically predicts a decou- 
pling of the gravitational potentials, ^ tp, such that 
the two potentials are independent functions of the four 
space-time coordinates, t and x: 



ip(r, x) = [1 + vj(t, x)] x 0(t, x) . 



(2) 



The model proposed in CCM attempts to link the depar- 
ture from GR with the growth of an effective dark energy 
relative to normal matter and radiation, whereby it was 
motivated that 



vo(t, x) = voqPbe(t, x)/p m (r, x), 
Pde(t) 



VOq- 



Pm(r) 



(3) 



The second line is obtained by expanding to first order 
in perturbations on the homogeneous background, indi- 
cated by the overbar. Because vo already multiplies first 
order perturbation quantities, the spatial dependence en- 
ters at second order and is neglected here. Truly, any 
theory of gravitational slip must be scale dependent in 
order to both satisfy solar system constraints at small 



length scales and allow for novel cosmological phenom- 
ena at extra-galactic length scales. For the purposes of 
our investigation, focusing on cosmology, vo is considered 
to be homogeneous. 

The background evolution is taken to be identical to 
ACDM, in which case 



— 



Wo cT ( 1 + z ) 



(4) 



Hence, the gravitational slip remains negligible until late 
times, when the cosmic acceleration becomes apparent. 
Our PPF cosmological model is described by the stan- 
dard set of ACDM parameters, plus voq. Informally, this 
model may be referred to as "hjACDM" . While there is 
no a priori reason that new gravitational physics should 
behave this way, this model provides a simple means to 
test for indications of new gravitational physics. If cos- 
mological data are found to be favoring a non-zero value 
for vo, an advanced theory of gravity beyond GR will be 
required to explain its value and rcdshift evolution. 

Our parameterization of the relationship between the 
gravitational potentials can be compared to modifica- 
tions of GR in the recent literature. Bertschinger & 
Zukin (2|| (hereafter BZ) adopt a notation with tp = $bz 
and <p = ^bz and parameterize 



*Bz(k, t) = 7 Bz(a)$Bz(k, t) 



(5) 



with 7bz(«) = 1 + (3a? where (3 and s > arc model 
parameters. In our language (|5|) takes the form <j) — 
7bz(o)-0, with 7bz(o) = 1/(1 + vo(a)) leading to 



vo(a) = -/3a s (l+ Pa 8 )- 1 



(6) 



In the case s = 3 and the limit of small /3, the 
paramctrization of BZ agrees with our model @ for 
(3 = —vooHa/^vi- As we discuss later, CMB anisotropy 
spectra obtained using the parameterization of equa- 
tion ^ in our code are fully consistent with the results 
of BZ. The differences in the CMB anisotropy spectra 
reported in BZ and CCM are due to a numerical error 
in the implementation of CCM, which we correct as dis- 
cussed below. 

Hu [28[ presents another parameterization with 



4' 



Hu 



Hu 



g(k,a) 



(7) 



with g(k,a) = g m (a)/(l + c 2 g k 2 H ) and g S n(a) = 

go \/ PdcAot / Ptot ^dc • In our notation the perturbations 
are now <j) = — $Hu and ip = $Hu- Ignoring the 
momentum dependence, which only dampens the post- 
Friedmannian modification below a certain scale, the cor- 
respondence is found to be 



vo(a) 



-2g(o) 
l + g(a) 



(8) 



Given the difference in rcdshift dependence between 
equations (01 and © for the g(a) used in Ref. [28j], a 
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simple relation between wq and go that is accurate at 
redshifts other than z = is not evident. Thus, a de- 
tailed comparison is not pursued. 

Jain & Zhang [25| present a post-Friedmannian pa- 
rameterization with 4> = r](k, t)tp. Ignoring the momen- 
tum dependence, rj is related to vj(o) through zu(a) = 
rj(a)^ 1 — 1. Since they do not present numerical results 
on the modification imposed by their parameterization to 
CMB anisotropics and other large-scale structure observ- 
ables, a comparison with their work is also not attempted 
here. 

In addition to modifications of gravitation that lead 
to <f> ^ if}, constraints on possible departures from the 
Newtonian inverse-square law at cosmological distance 
scales have also been considered [H, HH • It is not certain 
that such modifications satisfy the consistency relation 
for cosmological perturbations described in Ref. [ll[ ■ In 
any event, our study differs from such works in that no 
scale-dependent corrections to gravitation are imposed. 



III. IMPLEMENTATION OF PPF MODEL 

An implementation of the parametrized relationship 
between the potentials requires a consistent method of 
treating perturbations in the modified theory. To begin, 
the phenomenological relation between the gravitational 
potentials is defined in the conformal-Newtonian (longi- 
tudinal) gauge. In practice, however, we find it most 
convenient to pursue perturbation evolution in the syn- 
chronous gauge. Using the notation of Ref. [lj|, the met- 
ric perturbation variables in the two gauges are related 
as 



1 

'/ 



h - 
1 



f 6fj + H (h 
n(h + 6f) 



6fj 



(9) 
(10) 



where the dot indicates the derivative with respect to 
conformal time. In the standard GR case (w = 0), the 
perturbed Einstein equations, 



k ' 2 V - \nh = 47rGa 2 (5T 



h + 2Hh 



k 2 f] = 4irGa 2 (p + p)0 
2k 2 rj = -8nGa 2 STi, 



(11) 
(12) 
(13) 



are used to evolve the metric variables, where (p+p)9 = 
i&STj (for greater detail, see [13]). R is standard prac- 
tice to use the latter two equations (fT2l[T3T) for evolution, 
and apply the first equation (jll[) as a constraint. 

We presume a theory of modified gravitation in which 
the stress-energy tensor of matter and radiation is 
conserved, equation (|12p remains valid, but equations 
(fTTlfTB"]) are invalid. If the perturbed Einstein equations 
were assumed to remain valid, then the gravitational slip 
would necessarily imply the existence of new energy den- 
sity and pressure perturbations which are comoving with 



the baryonic and dark matter density perturbations. Fur- 
thermore, the non-zero gravitational slip introduces a 
modification to the perturbed, off-diagonal space-space 
Einstein equation 

a = -{2 + vj)Ha + {l + vj)i 1 -l2nGa 2 {p + p)(r/k 2 (14) 

where a = (h + 6r))/2fc 2 . The factor of 1/k 2 in the last 
term on the right hand side of (fT4")l corrects a typograph- 
ical error in CCM equation (8). Note that the shear due 
to matter and radiation is negligible at late times. This 
leaves a = — (2 + zu)Ti.a + (1 + ro)/? which is simply a 
restatement of (fTU)) in synchronous gauge. In order to 
maintain continuity from early times, when GR is valid, 
to late times, when the gravitational slip becomes signif- 
icant, equation (fT4Ij is used. In turn, 



q = v{l + m)-{2 + m)(Ha + Ha)+ib(Ti-Ha) 
-jL(l2irGa*(p+p)£). (15) 

Now, instead of evolving equations (fT2l[T3"|) with con- 
straint (fTTj) , equations (|12U14[) are evolved with the con- 
straint 



2k 2 c 



24irGa 2 



(p + p)6 



(16) 



derived from a and equation (|12p . Thus, the time evo- 
lution of 77, h and a, a are fully specified in a modified 
gravity model. 

The above set of equations, corresponding to the recipe 
Rl of CCM, has been implemented in a modified ver- 
sion of CMBfast [13] ■ As shown in Appendix [A] this 
prescription satisfies the consistency condition that was 
first derived in Ref. [ll| (equation 2 of BZ). This consis- 
tency suggests that our implementation, introduced by 
CCM, is not fundamentally different from that employed 
by BZ. Despite the reasons suggested in BZ for differ- 
ences between the two approaches, consistent results are 
obtained when the two post-Friedmannian parameteriza- 
tions are matched. 

Note that when presenting results here, the error in the 
numerical software used by CCM, involving a minus sign 
error in the w term of equation ()15p , has been corrected. 
Having rechecked the derivations and software, the re- 
vised results for the effect of w on the CMB anisotropy 
power spectrum are described below. 



IV. EFFECTS ON CMB ANISOTROPY 

The predicted CMB anisotropy spectra are compared 
with the WMAP 3-year data, and a likelihood function 
for the modified gravity variable vcq is constructed, using 
the likelihood code supplied by the WMAP team [3l|. To 
avoid duplicating a lengthy search of parameter space, 
the background cosmology is set to the 3-year WMAP 
maximum likelihood cosmology (hereafter WMAP3 ML: 
fl h = 0.0414, Q c = 0.196, fl A = 0.7626, h = 0.732, 
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FIG. 1: CMB anisotropy power spectra for different wo cos- 
mologies are shown with the binned WMAP 3-year data. All 
differences are localized to the low multipole moments: the 
spectra are normalized so that the higher multipole moments 
for all models are identical to the case of wo = 0, correspond- 
ing to the WMAP3 ML model. 
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FIG. 2: The predicted CMB temperature anisotropy 
quadrupole power is shown as a function of wo (solid curve). 
All other parameters are set to the WMAP3 ML cosmology. 
The central value and la upper bound of the WMAP 3 year 
data, 6C2/27T = 211 ± 860 |31(, are shown by the dashed and 
dot-dashed curves. 



r = 0.091, n s = 0.954), and for comparison to the 3- 
year WMAP plus SN Gold maximum likelihood cosmol- 
ogy (hereafter WMAP3+SNGold ML: fl h = 0.0454, fl c = 
0.2306, «a = 0.724, h = 0.701, r = 0.079, n s = 0,946) 
[34| . Figure Q] plots the multipole moments for several 
simulations with different values of zuq in the WMAP3 
ML cosmology. Normalizing the spectra to the observed 
amplitude of the acoustic peaks at t > 100, we see that 
only the lowest multipole moments, I < 30, are affected 
by wq. 

Furthermore, we observe that the large-angle 
anisotropy power grows as \tuq — 0.5 deviates from 



zero. For a closer look, the quadrupole moment is 
shown as a function of nj in Figure O This quadratic 
behavior can be understood as follows. The integrated 
Sachs-Wolfe effect contribution to the temperature 
anisotropy depends on <f> + ifi, which varies with w as 

<j> + i> = (3n7+ (r- l){2 + <w))H<t> (17) 

where T = 1 + d In cj>/d In a. Numerical results show that F 
decreases monotonically, becoming negative with increas- 
ing zuq- Upon squaring the above quantity, it is evident 
that \vj\ 3> 1 will lead to strong anisotropy in the low 
multipolcs, which are sourced at late times. Note that a 
similar trend, a local minimum of large-angle anisotropy 
power as a function of the post-Fricdmannian parameter, 
is seen in Hu's model, in Figure 1 of Ref. [2^]. BZ should 
see the same trendjff a wider range of (3 is searched than 
in Figure 5 of Ref. [26[ . 




FIG. 3: The likelihood of w due to the CMB, calculated 
using the November 2006 version of the WMAP likelihood 
code [3l[, is shown. All other cosmological parameters are set 
to the WMAP3 ML (solid curve) or the WMAP3+SNGold 
ML (dashed curve) model. The results include the TT and 
TE spectra. The results are consistent with wo — (ACDM), 
but favor positive values of wo- The locations of the primary 
(left) and secondary (right) likelihood peaks correspond to the 
range of wo for which the predicted quadrupole lies within the 
2a range of WMAP. 

The local minimum in the quadrupole amplitude helps 
explain the double-peaked structure of the WMAP like- 
lihood function, shown in Figure [3] The likelihood is 
suppressed at large values of \vj — 0.5 , due to the ex- 
cessive large angle anisotropy power. And while the low 
quadrupole moment reported by WMAP has been inter- 
preted to be indicative of new physics, in fact the WMAP 
likelihood docs not necessarily reward an anisotropy spec- 
trum that smoothly leads to low power at low £. The 
anisotropy spectrum with lowest quadrupole, at tuq = 
0.5, is suppressed relative to n7 = 0, 1.5 for the WMAP3 
ML parameters, as shown in Figure [3] In the case of 
the WMAP3+SNGold ML model, which has a slightly 
higher matter density, the primary peak in the likelihood 
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distribution is shifted towards a slightly higher value of 
vjq. To understand this behavior, note that in the case 
vo = 0, the strength of the ISW contribution depends on 
fi TO , with r w (f2 m [a])°- 55 . An increase in f2 m thereby 
diminishes the ISW contribution. Because the WMAP 
CMB data appears to prefer a moderate ISW contribu- 
tion, the low quadrupole notwithstanding, then the data 
should prefer a low matter density at wq = 0, as seen in 
Figure [3] 
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Figures [4] and [5] compare the results of BZ with our 
PPF implementation. To recreate BZ's results, a second 
CMBfast code was modified to evolve according to GR 
until z = 30. At that point, the parameter 



u 



3 1 



(18) 



is calculated, where w is the background equation of state 
(assuming that only matter and A contribute). Next, £ 
is set to remain constant for the rest of the calculation, 
whereupon the evolution is determined by the equations 




Hip 



(19) 
(20) 



and tp = 4>/-fBZ, 7bz = 1 + fla s , 77 = <p + Ha. Here, 
(JTnj) comes from equation (2) of BZ and (|2"D|) comes from 
the definition of a. This second code uses none of the 
Einstein equations. The results, shown in Figures |4] and 
OH are equivalent to our results produced using equation 
(|3]) if the appropriate relationship between w and 7bz 
is imposed. Appendix [X] demonstrates this equivalence 
analytically. 



FIG. 4: The equivalence of alternative techniques for evolv- 
ing cosmological perturbations under PPF gravitation is illus- 
trated by the effect on CMB anisotropy spectra. Lines show 
the spectra produced following our procedure (equations 1141 
I16[) with the parameterization ([4]) replaced by that proposed 
in BZ, zu(a) = — /3a s /(l + /3a") with s = 3. Symbols indicate 
the spectra produced following BZ's procedure (equations [18} 
I20|) . The agreement between the techniques is exact. 
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FIG. 5: The effect of the parameterization of gravitational slip 
on CMB anisotropy spectra is shown. Lines show the spectra 
produced using our parameterization, equation Q; symbols 
indicate spectra produced using the BZ's parameterization, 
4>/^ = 1 + Pa 3 with s = 3. To compare, the parameter f3 is 
set to /3 = — ■cuoQ.k/Qm- The agreement is excellent in the 
limit of small /3. 



V. THE GROWTH OF STRUCTURE 

The growth of density perturbations in baryonic and 
dark matter is affected by gravitational slip, as the po- 
tential produced by an overdensity no longer matches 
the potential responsible for gravitational acceleration. 
In GR, equations (fTTj) and (fL 7 !)) combine with the fluid 
conservation law (equation 29a in Ref. [l(|), 



5 = -a 




"I 



(21) 



to give the equation for the growth of non-rclativistic, 
prcssureless perturbations: 



6 + T-L5 = 4irGa 2 6p. 



(22) 



In the case w ^ 0, an independent differential equation 
for the matter density contrast is not readily available 
as equations f) 1 HI 1 3j> have been eliminated. CCM try to 
circumvent this obstacle, using (fT2l [T4")) and a, together 
with the assumption of negligible shear to arrive at 



S + HS = k 2 (l + m){Ha-Ti). 



(23) 



Unfortunately, that is as far as CCM can go without 
postulating another perturbed fluid SpnE to account for 
the theory's departure from pip . However, there is no 
barrier to evolving 8 — it is simply obtained numerically 
from ([2"T|) . 

The modified Boltzmann code CMBfast calculates 8 
as a function of time in the process of calculating the 
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FIG. 6: The matter density contrast as a function of scale 
factor is shown for different values of roo ■ The density contrast 
is normalized to a vjq = model. All models use WMAP3 
ML parameters. 



FIG. 7: The goodness of fit of the mapping of the growth rate 
between zuACDM and ACDM models, is illustrated. The 
percent difference in the density contrast for a model with 
zuo = 0.4, Q m = 0.3, and a ACDM model with tt m = 0.35. 
These cosmologies are predicted to be equivalent by equation 
(HU). 



anisotropy spectrum. Figure [6] illustrates the effect of 
zuq 7^ on the growth of long wavelength matter den- 
sity perturbations. This effect is monotonic with tuq > 
enhancing the growth at late times and wq < suppress- 
ing it. By comparing the rate of growth for models with 
different values of vjq and TO , an approximate mapping 
for the growth of 6 from njACDM to ACDM models has 
been obtained. Specifically, to a good approximation, the 
growth of structure in a wq ^ cosmology is equivalent 
to a ACDM cosmology with 

| ACDM = ^m|roACDM + 0.13g7 (24) 

In the parameter ranges —0.5 < vjq < 0.5 and 0.2 < 
^mlroACDM < 0.5, this fit is good to the 3% level for 
a > 0.2. (See Figure EJ). Note that m ^ docs not 
change the shape of the (linear) matter power spectrum. 

Figure [5] plots the predicted dimcnsionless linear 
growth rate, / = dhiS / dh\a, against a compilation of re- 
cent data. In the future, this may provide an alternative 
test of modified gravitation. Unfortunately, the current 
uncertainties in / are too large to provide a meaningful 
constraint on zuq. 



VI. WEAK LENSING 

Gravitational lensing phenomena depends directly on 
the sum of the two gravitational potentials, and is there- 
fore an excellent probe of gravitational slip. In the cos- 
mological setting, measurements of weak lensing of the 
pattern of galaxy clustering can be used to constrain vjq. 
Proceeding, the E-mode weak lensing convergence corre- 
lation function, £e f36| . is calculated, which requires the 




12 3 

z 

FIG. 8: The dimensionless linear growth rate / plotted as a 
function of z for different ujq cosmologies with the WMAP3 
ML background. The data point at z — 0.77 is due to 
Ref. [35|. The remaining, hybrid set of data points are taken 
from Table 1 of Ref. [j| , and are based on estimates of large- 
scale structure power spectrum growth at different redshifts. 

convergence power spectrum, given by [13, [H, [H, H(| 

p " (,, = i4 r^ p « h ( fc= ^") (25) 

»M = f*'*') 1 ^ (26) 

where x 1S the comoving radial distance with Xh is the 
horizon distance, r(x) is the comoving angular distance 
(x in the case of a flat universe), p(x) is the probability 
distribution of lensed sources, and refers to the present 
epoch. 

The lensing convergence depends on the power spec- 
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FIG. 9: The values of as resulting from the normalization to 
the WMAP 3-year anisotropy power spectrum, are shown as 
a function of zuo- All other cosmological parameters are set 
to match the WMAP 3 ML model. 
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FIG. 10: Model predictions of £e for different values of vjq. 
Data is taken from Table B.l of (4^. The background cos- 
mology is the WMAP3 ML. The power spectrum amplitude 
is determined by the normalization to WMAP. 



trum of metric perturbations (j> + tp. At linear scales, 
the power spectrum can be expressed through the power 
spectrum of density perturbations 6 by making use of the 
growth function Z? ra = 8(z)/S(z = 0) to account for the 
redshift growth of perturbations and the Poisson equa- 
tion. Using the relation between <p and ip, we can write 



P<f, + Tp(k, z) 



Ho 

" m ' v V ck 
xPss(k, z 



-n 2 

• si m,0 



DM 
a(z) 

0), 



w(z) 



(27) 



where Pss is the power spectrum of density perturbations 
today. Because the length scales probed by the weak 
lcnsing measurements extend into the non-linear regime, 
this description must be extended in order to compute 
the non-linear matter power spectrum. Since we do not 
have a complete description of non-linear evolution un- 
der modified gravity, here we employ the same technique 
as used for the usual GR predictions of lensing statis- 
tics. We use the fitting function from Peacock and Dodds 
[4lj | (PD) to get Pgg nhn from the linear power spectrum 
calculated with our modified CMBfast code. Note that 
PD's factor gpjj(Q, a) = 5(a) /a is meant to compare the 
growth history to that of an fl m = 1 model so 6(a) 
normalized at an early redshift is used. Specifically the 
density contrast is set to 6 (a = 0.01) = 0.01 for all wq 
models. The linear power spectrum amplitude is set by 
the normalization of the CMB anisotropy spectrum to the 
WMAP 3 year data. In Figurc[9]wc highlight the relation 
between as as a function of zuq and note that an increase 
in wq has a similar effect of enhancing weak lcnsing as 
does an increase in as in standard ACDM cosmologies. 
For wo = 0, the GR case, we recover the WMAP3 ML 
preferred value of cr 8 = 0.76 (34| . 

From the convergence spectrum, we can then find the 



modes of the correlation function |44J, |45|, |46l |4 

£f(0) 
£e(*) 



dkkP K (k)J {k6) 

-!- / dkkP K (k)J 4 {k9) 
27T Jo 

m±m ' (28) 



where Jiil) are Bcsscl functions of the first kind. Figure 
[TO] plots \ E (6) for different vj models with the WMAP3 
ML cosmology. Varying vjq results in a multiplicative 
shift in £e. 

Figure [IT] plots the likelihood derived from the weak 
lcnsing shear correlation function as measured by the 
Canada- France-Hawaii Telescope Legacy Survey [481 ] and 
presented in Table B.l of Fu et al. [43|]. The source dis- 
tribution p(x) is taken from equation (14) and Table 1 of 
Ref. [43j]. Unlike the CMB anisotropy spectra, the weak 
lensing shear correlations functions do not "bounce" with 
ujq, so the resulting likelihood function, shown in Fig- 
ure [11] is narrower and peaks at wq ~ 0.3. This result 
excludes vj§ = 0, but is roughly consistent with the pri- 
mary peak in the CMB likelihood, shown in Figure [3] 

The exclusion of wq = based on CMB and weak 
lensing data in the WMAP3 ML model is not signifi- 
cant at this stage. Only a single parameter in a multi- 
dimensional parameter space has been searched. Indeed, 
a positive value of wq is expected since we have nor- 
malized the underlying density power spectrum to the 
WMAP 3 preferred value of as, whereas analyses of the 
weak lensing data have found a preferred normalization 
value with as ~ 0.8, which is slightly higher than the 
CMB value. The disagreement between these two val- 
ues is now resolved by a non-zero value for wq. While 
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FIG. 11: The likelihood of different vjq models according to 
the £e data collected by the Canada-France-Hawaii Telescope 
Legacy Survey and presented in Table B.l of Ref. [3] is shown 
with parameters set by the WMAP3 ML model (solid line) as 
well as the WMAP3+SNGold ML model (dashed line). 



such a difference between CMB and low-rcdshift matter 
perturbations is generally the very type of behavior we 
should look for in modified gravitation scenarios, the fact 
that we do not perform a joint combined analysis of CMB 
and weak lensing data precludes us from making a strong 
statement. 

We expect a degeneracy in the effect of wq and £l m h 2 
on the weak lensing predictions, as both parameters con- 
trol the growth rate of fluctuations. Specifically, increas- 
ing wq enhances fluctuation growth, as seen in Figure [6l 
just as does increasing the abundance of clustering mat- 
ter. By raising fl m h 2 the likelihood is expected to peak 
at lower values of vdq- This is precisely illustrated in 
Figure 1111 wherein the weak lensing likelihood for wq is 
shown for cosmological models with parameters set by 
the WMAP3+SNGold ML model, which has a slightly 
higher matter density. This new likelihood is indeed con- 
sistent with lower values of hjq, including zero. 



VII. CMB AND LSS CROSS-CORRELATIONS 

A deviation from GR leaves an imprint on the cross- 
correlation between the CMB and large-scale structure. 
In order to identify such a signal, we define the two-point 
angular cross-correlation between the temperature ISW 
anisotropy and the dark matter fluctuation as [ill [5(3, [EH, 

m 

C x (6) = (Ajswd^Lss^)), (29) 

where the angular brackets denote the average over the 
ensemble and 9 = |t\ — 72 1. For computational purposes 
it is convenient to decompose C x (9) into a Legendre se- 



ries such that, 

00 0/ , 1 

C X (e)=^2^C^P l (co S (9), (30) 

1=2 

where Pi (cos 9) are the Legendre polynomials an d C, x is 
the cross-correlation power spectrum given by [5(3, HJ, [H, 

m 

C x =A^JfAlir w (k)lt ss (k), (31) 

where is the primordial power spectrum. The inte- 
grand functions l[ (k) and // jSS (fc) are defined respec- 
tively as 

Il LSS (k) = b J $(z)5 k (z)ji[kr(z)}dz, (33) 

where 4>k and 8 k are the Fourier components of the gravi- 
tational potential and matter perturbation, respectively; 
$ is the galaxy survey selection function; ji [kr(z)] are the 
spherical Bessel functions; r(z) is the comoving distance 
at redshift z and k(t) = k(t)cIt is the total optical 
depth at time r. 

A change in w could therefore change not only the 
value of l( sw but also its sign with respect to /; iS ' s . 
Direct measurements of the cross-power spectrum C ; 
are more robust for likelihood parameter estimation since 
these data would be less correlated than measurements 
of C x (9). Therefore, the cross-power spectrum C x is 
computed for different values of wq assuming a galaxy 
survey with a selection function as 

$(z) - z 2 cxp[-(z/z) 1 ' 5 ] (34) 



0.12, 




10 10 
multipole moment (1) 



FIG. 12: The cross-correlation angular power spectrum be- 
tween CMB temperature (ISW) and galaxy distribution is 
shown as a function of wo- A value of 0.5 < zuq < 1.5 changes 
the sign of the ISW, resulting in a negative cross-correlation, 
in conflict with observations. 
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O Gaztanaga et al. 



□ Cabre et al. 
U.U1 . a Giannantonio et al. 

"iii i , i , i i i : 

0.1 1 
z 

FIG. 13: The cross-correlation amplitude at = 6 degrees is 
shown as a function of redshift for different samples of LSS 
and WMAP data. The curves show the expected correlation 
with vjo varied. The data plotted here comes from the compi- 
lation in Ref. [6(j, making also use of data from Gaztanaga et 
al. [6ll (c ircles). Cabre et al. [f^] (squares), and Giannantonio 
et al. [60| (triangle). 



where z, the median redshift of the survey, is z = 0.25. In 
recent years, the WMAP temperature anisotropy maps 
have been cross-correlated with several surveys of Large 
Scale Structure (LSS) distributions and a positive corre- 
lation signal has been detected [H, |H, H3, HI, [H, [6(| • 
As seen in Figure [T^l values of wq > 0.5 would result 
in anti-correlation and a negative angular spectrum, in 
disagreement with current observations at the ~ 2er level. 

In order to obtain a more quantitative result, the ISW- 
Galaxy theoretical correlation function for different val- 
ues of w is compared with the current data as reported in 
Ref. [H under the assumption of the WMAP ML model. 



The data and some model curves under modified grav- 
ity are shown in Figure 1131 The likelihood is reported 
in Figure [Ml consistent with no indications for modified 
gravity. This constraint should be considered with cau- 
tion, since the full parameter space has not been explored 
and the ISW signal may certainly be highly sensitive to 
changes in parameters such as the matter density and the 
Hubble constant. 

Qualitatively, the ISW-Galaxy theoretical correlation 
may be understood by considering the path of a CMB 
photon through a region on the sky containing an abun- 
dance of galaxies. As the photon passes through the col- 
lective gravitational potential of the clustered galaxies, 
the decay of the gravitational potential in a ACDM cos- 
mology with T37o = leads to the blueshifting of the pho- 
ton. A negative value of w§ suppresses the clustering and 
causes the photon to be further blueshifted as it climbs 
out of the weakening gravitational potential. The coin- 
cidence of the hot spot in the CMB resulting from the 
blueshifted photons, with the abundance of galaxies on 
the sky explains a positive cross-correlation. Likewise, a 
large wq enhances clustering and reduces the blueshifting 
as the photon climbs out of the growing gravitational po- 
tential. In extreme cases, the enhanced clustering causes 
the photon to be redshifted, thereby leading to a negative 
cross-correlation. The degeneracy between vjq and fl m h 2 
suggests that a preference for vjq < seen in Figure [HI 
will be enhanced by increasing fl m h 2 : the slower decay 
of the gravitational potential (T — 1 in equation 1 17p due 
to the increased matter density is compensated by the 
reduced clustering with larger zuq. Figure [Til also shows 
the fit to current ISW-galaxy cross correlation data using 
the WMAP3+SNGold ML model parameters, for which 
the matter density is slightly higher than the WMAP3 
ML model. Indeed, the peak of the likelihood shifts to 
even more negative values, as expected. 




FIG. 14: The likelihood of different ujq models using current 
ISW-galaxy cross correlation data is shown. All other cosmo- 
logical parameters are fixed to the WMAP3 ML model (solid 
line) or the WMAP3+SNGold ML model (dashed line). 



VIII. SYNTHESIS 

We have examined the consequences of the proposed 
PPF model of gravitational slip for predictions of the 
cosmic microwave background anisotropy, the growth 
of large scale structure, weak lensing, and the ISW- 
galaxy cross-correlation. The likelihood distributions in 
wq using the WMAP3 ML model parameters are sum- 
marized in Figures [TBI and using the WMAP3+SNGold 
ML model parameters in Figure [111 The weak lensing is 
clearly the most sensitive to wq : followed by the ISW- 
galaxy cross-correlation, and then the CMB. 

In the case of the WMAP3 ML model parameters, 
shown in Figure 1151 the overlap between the three dis- 
tribution functions is strongest near wq ~ 0.2. However, 
there is no concordance due to the tension between the 
values vjq > indicated by the CMB and weak lensing 
data, and the values centered on zero for the ISW-galaxy 
cross-correlation. This tension appears to be related to 
the disagreement in the best-fit erg values derived from 
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measurements of the CMB and large scale structure. In 
fact, the additional degree of freedom introduced by Wq 
seems to be justified by the sharp peak in the weak lcns- 
ing distribution. 

To find concordance among the three observational 
constraints, we increased the matter density. As we ar- 
gued, this should have the effect of shifting the weak 
lcnsing distribution towards lower values of vjq, with 
the primary concern to determine whether vjq — is 
allowed. Increasing the matter density also has the ef- 
fect of moving the CMB distribution towards higher val- 
ues, away from zuq = 0, and moves the ISW-galaxy 
cross-correlation towards more negative values, also away 
from n7 = 0. Using the WMAP3+SNGold ML model 
parameters, as shown in Figure 1161 the overlap is im- 
proved. The tension between weak lensing and ISW is 
somewhat relaxed, with a joint likelihood that allows 
vjq = 0. Whether wq = is preferred remains to be 
demonstrated. A full, multi-dimensional parameter space 
analysis is planned for a subsequent paper. 



8 - m WMAP3 + SNGold ML 

1 1 
i i 
i i 




FIG. 16: The likelihood distributions for wo, based on the 
CMB (solid curve), weak lensing (dashed), and ISW-galaxy 
cross correlation (dot-dashed) using the WMAP3+SNGold 
ML model parameters is shown. 
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FIG. 15: The likelihood distributions for wo, based on the 
CMB (solid curve), weak lensing (dashed), and ISW-galaxy 
cross correlation (dot-dashed) using the WMAP3 ML model 
parameters is shown. 



IX. DISCUSSION 

A PPF formalism for modified gravity has been pro- 
posed to describe possible departures from cosmological 
predictions under GR. The model consists of a back- 
ground cosmology, a parameterized relationship between 
the gravitational potentials, and an implementation of 
new evolution equations in the absence of the perturbed 
Einstein's equations. For simplicity we have focused on 
the case in which the background cosmology evolves as 
a standard, ACDM universe. The gravitational slip is 
chosen to evolve in proportion to the dominance of an 
effective dark energy density over matter. The constant 
of proportionality introduces a new parameter, wq. The 



new evolution equations are summarized by equations 

We have made a comparison between our proposed 
PPF model of modified gravitation and other models 
that have recently appeared in the literature. Our im- 
plementation satisfies the consistency relation proposed 
by Bertschinger [111 ]. Furthermore, the equivalence with 
the model of BZ has been demonstrated, after accounting 
for the correspondence between w{z) and 7bz(o), using 
analytical calculations and numerical results. 

We remark that many different parameterizations for 
the departure from GR on cosmological scales have been 
introduced in the literature. This variety may be use- 
ful, as no single model of non-Einstein gravitation stands 
out. However, it would also be useful if the data are 
compared within a single parameterization, similar to the 
post-Newtonian parameter 7ppn that is used to measure 
the amount of spacctime curvature per unit mass in tests 
of gravitation in the solar system. In that regard, our 
parameter corresponds to vo s=s 1 — "fppN in the limit of 
weak departures from GR. When w is positive, then the 
gravitational potential determined by geodesic motion is 
greater than the potential inferred from the distribution 
of matter via the Poisson equation. 

An alternative view of the gravitational slip is that it 
effectively introduces a new source for the perturbed off- 
diagonal space-space Einstein equation. If we infer the 
existence of some new density perturbations, defined so 
as to balance the perturbed time-time Einstein equation, 
then the new source takes the simple form iro^ptot- 

We have corrected an error in CCM, and proceeded 
to examine the cosmological consequences of zuo =/= 0. 
We hope our study clarifies that: (a) there is a tension in 
current data that may be suggestive of PPF modifications 
of GR; (b) different parameterizations produce the same 
results at the end as long as certain consistency relations 
are satisfied; and (c) it would be useful to establish one 
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or two (and same) parameters from future data to finely 
test the departures from GR at cosmological scales. 



APPENDIX A: EQUIVALENCE BETWEEN 
MODELS OF PPF GRAVITATIONAL SLIP 

Bertschinger and Zukin [26| propose that 4> evolves ac- 
cording to equation (fl9|) after z = 30, with constant cur- 
vature perturbations, Q. To show that their model is con- 
sistent with roACDM, their evolution equation is recast 
as 



Hip 



V 



H(l-H) 
2 fj 
3H(l + w) 



(Al) 



where the translation from conformal-Ncwtonian (longi- 
tudinal) gauge to synchronous gauge, 



4> = rj — Ha 
ip = a + Ha , 



(A2) 
(A3) 



is used in the second equality. Also, in a matter- and 
A-filled universe, H(l ~ H/H 2 ) = §«(1 + w) where w is 
the background equation of state. 

For consistency with the roACDM model, the evolution 
of rj, fj according to equations (|14j|16p must agree that ( 
as given in equation (|A1|) is a constant. Hence, assuming 
the zero-i perturbed Einstein equation (equation 21b of 
[3) then 



= V 



2 

H9 



4iTGa 2 {p + p)9 
k 2 H{l+w) 



where the background Fricdmann equation is used to get 
the second equality. Taking a derivative with respect to 
conformal time, 



C = V + 



He ne 



(A4) 



Again using the zero-i Einstein equation, and the equa- 
tion for 9 derived from stress-energy conservation (equa- 
tion 29 of [l3]), 



= -H(l-3w)9- 



1 



t-e+SElSPtfs-h**, (as) 

w 1 + w 



then the time derivative of the curvature perturbation is 

, 5p S 



( = H( 



Spl 



(A6) 



Since the dominant perturbations at late times {z < 30) 
are due to baryonic and dark matter that have negligible 
pressure and shear perturbations, then £ w 0. Therefore, 
tuACDM is consistent with the evolution equation (jAlj) 
with constant £. 
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